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1. Introduction 

One of the most exciting areas in applied mathematics is the study of the dynamics 
associated with the propagation of information. Coherent structures hke solitons, 
kinks, vortices, etc., play a central role, as carriers of energy, in many nonlinear 
physical systems [l7|. Solitons represent a rare example of a (relatively) recently 
arisen mathematical object which has found successful high-technology applica- 
tions jlj]. The nature of the system dictates that the relevant and important effects 
occur along one axial direction. Interplay between nonlinearity and periodicity is 
the focus of recent studies in different branches of modern applied mathematics 
and nonlinear physics. Applications range from nonlinear optics, in the dynamics 
of guided waves in inhomogeneous optical structures and photonic crystal lattices, 
to atomic physics, in the dynamics of Bose-Einstein condensate (BEC) droplets 
in periodic potentials, and from condensed matter, in Josephson-junction ladders, 
to biophysics, in various models of the DNA double strand. Analysis and model- 
ing of these physical situations are based on nonlinear evolution equations derived 
from underlying physics equations, such as nonlinear Maxwell equations with peri- 
odic coefficients [37]. In particular, the systems of 2nd-order NLS equations, both 
continuous and discrete, were applied in nonlinear physics to study a number of ex- 
perimental and theoretical problems. Spatial non-locality of the nonlinear response 
is also naturally present in the description of BECs where it represents the finite 
range of the bosonic interaction. Demands on the mathematics for techniques to 
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analyze these models may best be served by developing methods tailored to de- 
termining the local behavior of solutions near these structures. The discreteness of 
space i.e., the existence of an underlying spatial lattice is crucial to the structural 
stability of these spatially localized nonlinear excitations. 

During the early years, studies of intrinsic localized modes were mostly of a 
mathematical nature, but the ideas of localized modes soon spread to theoretical 
models of many different physical systems, and the discrete breather concept has 
been recently applied to experiments in several different physics subdisciplines. 
Most nonlinear lattice systems are not integrable even if the partial differential 
equation (PDE) model in the continuum limit is. While for many years spatially 
continuous nonlinear PDE's and their localized solutions have received a great 
deal of attention, there has been increasing interest in spatially discrete nonlinear 
systems. Namely, the dynamical properties of nonlinear systems based on the inter- 
play between discreteness, nonlinearity and dispersion (or diffraction) can find wide 
applications in various physical, biological and technological problems. Examples 



are coupled optical fibres (self-trapping of light) [l, J^, [l^, 26], arrays of coupled 
Josephson junctions [s^, nonlinear charge and excitation transport in biological 
macromolecules, charge transport in organic semiconductors (4ol |. 

Prototype models for such nonlinear lattices take the form of various nonlinear 
lattices [6|, a particularly important class of solutions of which are so called discrete 
breathers which are homoclinic in space and oscillatory in time. Other questions 
involve the existence and propagation of topological defects or kinks which math- 
ematically are heteroclinic connections between a ground and an excited steady 
state. Prototype models here are discrete version of sine-Gordon equations, also 
known as Frenkel-Kontorova (FK) models, e.g. [3]. There are many outstanding 
issues for such systems relating to the global existence and dynamics of localized 
modes for general nonlinearities, away from either continuum or anti-continuum 
limits. 

In the main part of the previous studies of the discrete NLS models the dis- 
persive interaction was assumed to be short-ranged and a nearest-neighbor ap- 
proximation was used. However, there exist physical situations that definitely can 
not be described in the framework of this approximation. The DNA molecule con- 
tains charged groups, with long-range Coulomb interaction 1/r between them. The 
excitation transfer in molecular crystals and the vibron energy transport in 
biopolymers [IH] are due to transition dipole-dipole interaction with depen- 
dence on the distance, r. The nonlocal (long-range) dispersive interaction in these 
systems provides the existence of additional length-scale: the radius of the disper- 
sive interaction. We will show that it leads to the bifurcating properties of the 
system due to both the competition between nonlinearity and dispersion, and the 
interplay of long-range interactions and lattice discreteness. 

In some approximation the equation of motion is the nonlocal discrete NLS 



lUr, 



+ U- 



n G Z , 



(1) 



where the long-range dispersive coupling is taken to be either exponentially Jn = 
jg-/3|"l with /3 > 0, or algebraically J„ = J|n|~* with s > 0, decreasing with the 
distance n between lattice sites. In both cases the constant J is normalized such 
that X^^i Jn = ^1 for all /? or s. The parameters (3 and s are introduced to cover 
different physical situations from the nearest-neighbor approximation (/? — > oo, s — > 
oo) to the quadrupole-quadrupole (s = 5) and dipole-dipole (s = 3) interactions. 
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The Hamiltonian H and the number of excitations N 

= ^ ^ Jn-m\Un " Um]'^ ~ ^ '""l"^' ^ = X] (2) 

are conserved quantities corresponding to the set of ([T]). 

It should be also noted that the derivation of a discrete equation from the Gross- 
Pitaevskii equation produces at the intermediate step a fully nonlocal discrete NLS 
equation for the coefficients of the wave function expansion over the complete set 
of the Wannier functions. Further reduction to the case of the only band with 
the strong localization of the Wannier functions (the tight-binding approximation) 
leads to the standard local DNLS equation. Recently Abdullaev et al. [2] extended 
this approach to the case of periodic nonlinearities and derived a number of non- 
integrable lattices with different nearest-neighbor nonlinearities. 

In this paper, we study the discrete nonlinear Schrodinger equations on the lattice 
Z (DNLS) with nonlocal interactions of forms 

lUn = '^aj/\jUn + f{\Un\^)Un, nGZ (3) 

where ti„ G C, AjU^ := Un+j + Un-j — 2un are 1-dimensional discrete Laplacians 
and it holds 

(HI) / G C(IR+, IR) for R+ := [0, oo), /(O) = and aj G R with Wj\ < oo. Moreover, 
there are constants s > 0, ^ > 1, ci > 0, C2 > and f > such that 

\f{w)\ <ci{w' + 1), C2{w''+^ -I) < F{w), nF{w) -r < f{w)w 

for any w >0, where F{w) = J f{z)dz. Furthermore, limsup^^g fiw)/^^ < oo 

_ ^ 

for a constant s > 0. 

Of course we suppose that not all aj are zero. Note any polynomial f{w) = piw + 
■ ■ ■ +psw'^, s G iNi with ps > satisfies (HI). Furthermore, ^ can be rewritten into 
a standard form 

lUn = ^ ai^_n\{Um - Un) + f{\Un\'^)Un, n G Z. (4) 

It is well known that ^ conserves two dynamical invariants 

|iinP — the norm, 



E 

nez 



9 X] "I 

m—n\ 

\Um-U„\ +F{\Ur^ 



2 



the energy. 



Differential equations with nonlocal interactions on lattices have been studied 
in SJlJi pLMJll, Q, m, [sol], while DNLS (discrete nonlinear Schrodinger) in 
10, U, 3, 22; 28l|. Nowadays it is clear that a large number of important models 
of various fields of physics are based on DNLS type equations with several forms 
of polynomial nonlinearities starting with the simplest self-focusing cubic (Kerr) 
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nonlinearity, then following with the cubic onsite nonlinearity relevant for Bose- 
Einstein condensates, then with more general discrete cubic nonlinearity in Salerno 
model up to cubic-quintic ones (see Ul] for more references). 

We are interested in the existence of traveling wave solutions Un{t) = U{n — ut) 
of ([3]) with a quasi periodic function U{z), z = n — vt and some v ^ Q. 

First, we introduce a function 



:= — aj siv? 



Remark 1 : Clearly ^> G C(IR \ {0}, R), $ is odd, = for any A; G Z \ {0}, 

and ^ as |x| ^ oo. If ^ j\aj\ < oo then <I> e C{R, R) and if J^kil < ^ 

then ^ E C'^C'R, IR)- Consequently the range TZ^ := $(IR \ {0}) is either an interval 
[—R,R] or {—R,R) here with possibility R = oo (see Section [2.41 for concrete 
examples). 

Now we can state the following existence result. 

Theorem 1.1 : Let (HI) hold and T > 0. Then for almost each u £ R \ {0} 
and any rational r G Q n (0, 1), there is a nonzero periodic traveling wave solution 
Un{t) = U{n - vt) of ^ with U G C^(IR, C) and such that 

U{z + T)= e^^"U{z), Vz e R . (5) 

Moreover, for any G IR \ {0} there is at most a finite number of fi,f2, ■ ■ ■ ,fm G 
(0, 1) such that equation 

-u = ^(^^{rj + k)^ 

has a solution k G Z. Then for any r G (0, 1) \ {fi, f2, • • • , f^} there is a nonzero 
quasi periodic traveling wave solution Unit) = U{n — vt) with the above properties. 
In particular, for any \v\> R and r G (0, 1), there is such a nonzero quasi periodic 
traveling wave solution. 

When a nonresonance condition of Theorem 11.11 fails . then we have the following 
bifurcation results. 

Theorem 1.2: Suppose f G C'^{R+,R) with /(O) = 0. // there are n G (0,1), 
I' G TZ^ \ {0} and T > such that all solutions ki,k2, - ■ ■ , fcmi G Z of equation 

-u = ^(^^{n+k)^ 

are either nonnegative or negative, and mi > 0. Then for any e > small there are 
mi branches of nonzero quasi periodic traveling wave solutions Un,j,e{t) = Uj,s{n — 
Vet) of ^ with Uj^e e C^(R, C), j = 1, 2, • • • , mi, and nonzero velocity satisfying 
Uj^e^z + T) = e^'^''"^*C/j-.e(z), Vz G IR along with ^ v and Uj^e ^ uniformly on 
IR as e ^ 0. 

Remark 2 : If Oj > for all j G IH, then the assumptions of Theorem 11.21 are 
satisfies for any i/ G TZ^ \ {0} such that ^^~^{—i^) \ ^ / 0, and so there are 
bifurcations of quasi periodic traveling waves in the generic resonant cases. On the 
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other hand, Hue TZ^ \ {0} with 
for any r G (0, 1). 



-v) C J. then Theorem II. II is apphcable 



Theorem 11.21 is a Lyapunov center theorem for travehng wave solutions. Similar 
results are derived in 23| for Fermi-Pasta-Ulam lattices. 

We also discuss in Section H] the extension of these results of ([3]) on the lattices 
and 0, [HI, 112, HI] . The final Section [5] is devoted to travelina; wave solutions 
of more general forms than above ^] . 



2. Existence of Traveling Wave Solutions 

In this section, we study the existence of traveling wave solutions of the form 
Un{t) = U{n — i/t), i.e. we are interested in the equation 

- mU\z) = ^ ajd.Uiz) + f{\U{z)\^)Uiz) , (6) 

where z = n—vt, / and djU{z) := U{z+j) + U{z—j) — 2U{z). We are interested 
in the existence of quasi periodic solutions U (z) of ([6]) stated in Theorem 11.11 



2.1. Preliminaries 

In this subsection we recall some results from critical point theory of [27]. Let H 
be a Hilbert space and let J G C^{H, IR). Suppose H = Hi (B H2 for closed linear 
subspaces, and let 61,62,- •• be the orthonormal basis of Hi. Let us put H^ := 
spanjei, 62, • • • , 6^} and i/„ := H^ ® H2. Let Pn be the orthogonal projection of 
H onto Hn- Set Jn := J/Hn - the restriction of functional J on subspace Hn - and 

so VJnix) = PnVJix) if X G F„. 

Definition 2.1: If there are two positive constants a and /3 such that 

Jix) > 
J{x) > a 
J{x) < 
J{x) < —a 

then J is said to satisfy the local linking condition at 0. 

Definition 2.2: We shall say that J satisfies the Palais-Smale (PS)*-condition if 
any sequence {x„}neiH in H such that Xn G Hn, J{xn) < c < 00 and PnVJ{xn) = 
^Jn{xn) as n ^ 00, possesses a convergent subsequence. 

Now we can state the following theorem of [27I] which we apply. 

Theorem 2.3 : Suppose 

ih) JG C^(i?, IR) satisfies (PS)* -condition. 
{I2) J satisfies the local linking condition at 0. 
(/s) Vn, Jn{x) — > —00 as \\x\\ 00 and x G Hn- 

{I4) V J = A + C for a bounded linear self-adjoint operator A such that AHn C Hn, 
Vn G iNi and C is a compact mapping. 

Then J possesses a critical point x with \J{x)\ > a. 



\/x £ {x £ Hi \ \\x\\ < (3} , 
Vx G {x G -fTi I ||x|| = P} , 
Vx G {x G F2 I \\x\\ < /?} , 
Vx G {x G -^2 I Ikll = /?} , 
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Remark 1 : If is an indefinite nondegenerate critical point of J, then J satisfies 
the local linking condition at 0. 



2.2. Proof of Theorem [7711 

In this section, we use Theorem 12.31 to prove Theorem ll.il Without loss of gener- 
ality, we set T = 27r. We suppose v > 0, the case u < can be handled similarly. 
First, we identify C with IR^ in this section. Let r £ (0, 1) be fixed. Next, we consider 
real Banach spaces 

:= {[/ e lL(R, C) I U{z + 2tt) = e^^"U{z), Vz G r} 

for s > 1. Clearly [/ G if and only if U{z) = e'^^^V{z) for some V G L'^ := 
C). Consequently Ui{z + ci)U2{z + C2) is 27r-periodic for any ci, C2 G IR and 
Ui,U2 G hence is 27r-periodic. So we consider the norm on Lf. like on 

L^. In particular, we have 

y G l2 4^ V{z) = Y,Vk e^'+^>\ VkeCY, iVkf < 00 . 
feez fcez 

Let 

Xr := Wy^^\S^\C) = lvGLl\ V{z) = ^yfce(^+^)^\ \VkW + k\ < oo \ , 

I feez fcez J 

Yr := W,'^\S^\C) = \veLl\ V{z) = ^ 14 e(^+'=)^\ ^ |yfc|2(r + A:)^ < ooj . 

I fcez fcez J 



Note r + /c 7^ for any k £ Z. Clearly Yr C Xr C Lj. We consider Lf., Xr and Yr 
as real Hilbert spaces with inner products 

{V, W) := 2tt^ V = K / V{z)W{z)dz , 

fcez 

(y, W)x^ := 2^3f? FfcWfclr + k\ 
fcez 

(y, W)y„ := 2tt^ Y1 ^kW~k{r + kf 

fcez 



for V{z) = Vke^''+''>' and W{z) = ^ Wfc 6^'"+'=)^*. 

Clearly ||[/||l2 = ||[/||i2 < ri||[/||x., VC/ G X,. and < = pil^^^' 

VC/ G Yr for ri := min vT^-^}- The following result is well-known 271. l34[. 

Lemma 2.4: For each J>l,Xr is compactly embedded into L,^. 

On Xr, we consider a continuous symmetric bilinear form 



fcez 
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Note, if [/ G and F G Tr, then 

2tv 



/ iU{z)V{z)'dz = Br{U,V). 
Jo 



Now we consider a real functional 



dz 



= U) + |g a,-|C/(z + j) - ^(z)|2 - F(|[7(z)p)| 

on Xj.- Then 1^ G C^{Xr, R) and for U G X^, ^ G 1^, we derive 

DIr{U)V = 

23? <( ^ Li;7(z)y(I) ' - \^ajdjU{z) + f{\U{z)\^)U{z)\v{£)\ dz 



If [/ G Xr is a critical point of then 



2-K 



mU{z)V{z)'- I Y,a,djU{z) + f{\U{z)\'')U{z) \ V{z) \dz)=Q 



for any V ^Y^. Replacing V with «y in ([7]), we obtain 



(7) 



fj i^iU{z)V{z)' - i^jijdjU{z) + f{\U{z)\'')U{z)^ V{z)^ dz = 

for any V G Yj.. This means that [/ is a weak solution of ([6|). Then a standard 
regularity method shows 3j] that {/ is a C^-smooth solution of ([6]). 
Now we split = X_|- © X_ for 



-1 



X_ := <^y(z) = E ^fce 



(r+fc)zi 



X+:=\v{z) = J2Vke^'' 



k=—oo 



k=0 



Clearly if [/ = C/+ + C/_ then fi,(C/, [/) = 2 - ||C/-||xJ- 

Next, let us define Kr ■ Lf. Xj. as 



2tt 



{KrH, V)x^ ■■= 23f? / H{z)V{z)dz, W € Xr . 
Jo 



Then 



KrH = E 



2H, 



r + k\ 



fc „{r+fe)iz 



(8) 



September 10, 2009 17:14 Applicable Analysis revisedlatticeFeckanRothos 



Michal Feckan and Vassilis M. Rothos 



and so Kj. is compact. To study VIr{u), we introduce the mapping '• -^r — ^ -^r 
defined by 



2tt 



i^riU), V)x^ ■■= / fi\Uiz)\')Uiz)V{z)dz 
Jo 

for any V £ X.^. By Lemma [2. 4| the Nemytskij operator U f {\U {z)\'^)U {z) from 
Xr to L1 is continuous. Using ([8]), we get 

^r{U)=Krf{\U\^)U . 

Hence '■ Xr X^ is compact and continuous. 
Lemma 2.5: Under (HI) it hods D^r{G) = 0. 
Proof : There is a constant C3 such that 

I/HI <c^{w + w') 

for any to > 0. Then by Lemma 12.41 '^^ derive 

\fm^)U\li = / f{\U{z)\'f\Uiz)\'dz 



< 2c 



JO 

■I I {mz)f + |C/(z)| W)) dz < d (||[/||1 + WUf^^'f 
for a constant C4 > 0. Hence 

K^,(c/),F)xJ < 2\\fi\u\^)u\\LfMv\\Li < C5 (||c/|ll + ||[/|||+i) ||y|k. 

for a constant C5 > 0. This imphes 

||*r(?7)||x. < C5 (||?7|||^ + ||t/||^^+l) ,yUeXr. 

Since ^'r-(O) = and s > 0, we get D'^r{0) = 0. The proof is finished. □ 
Finahy, define Cr Li ^ L^. as 

CrU := ''^ajdjU{z) . 

Then 

VIr{U) = {2vl+ - 2ul_ - KrLr - (U) (9) 

for the identities I± : X± X±. Clearly 

Ar := 21^1+ - 2vl^ - KrCr 
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is a self-adjoint bounded operator Ar : Xr Xj. satisfying 



ArU = 2^ 1/ sgn (r + A:) + 



I — > aj sm 

\r + k\ ■' 



r + k 



Uke 



Consequently, the spectrum a{Ar) of Ar is given by 



a{Ar) = {2 sgn (r + fc) (i/ + ^>(r + A;)) | A: G Z} . 



By Lemma l2.5| we get that under the assumption 



(10) 



is an indefinite nondegenerate critical point of If - V/r(0) = and Hess /^(O) = Aj. 
with ^ 17(74^) and Xr = Xi^r®X2,r with a{Ar/Xi^r) C (0, 00) and a{Ar/X2,r) C 
(—00, 0) where Xi^r, -^2,r are suitable closed linear subspaces of Xr- Note Xi^^ and 
X2^r are infinite dimensional, since $(r + fc) — s- as |fc| — > 00. Consequently by 
Remark [H under (|lUp , satisfies the local linking condition at in the sense of 
Definition 12. H i.e. condition (I2) of Theorem 12.31 is verified. 
We consider an equivalent scalar product (•, •)r on Xr such that 



{ArU,U)r = \\Uifr " ll^2||', Ui G U2 G X2, 



Note there is a linear isomorphism Kr : ^ -'^r such that 



{U, V)x^ = {KrU, V)r, \/U,\/V G Xr . 



Clearly Kr is self-adjoint and positive definite. Then 

IriU) = ^\\Ult - ^\\U2fr - F{\U{z)\')dz, 
VIr{U) = Vli - UI2 -Kr-^r, 
{VIr{U),V)r = DIr{U)V = v\\Vi\\l " i^ll^2||r 

-2K r f{\U{z)\^)U{z)V{z)dz. 
Jo 

Let X2^r = span {ei, 62, • • • } and are eigenvectors of Ar. Then we take Xn = 
span {ei, 62, • • • ,6„} © Xi^r for n > 3. So clearly AXn C Xn, i.e. condition 
of Theorem 12.31 is verified. Let Pn '■ Xr Xn be the orthogonal projection with 
respect (•, •)r- 

We suppose there is a sequence {Um}m(^M C Xr, Um G Xm and a constant c such 
that 



IriUrr,) < C and PmV Ir{Um) ^ . 
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Then for m large we get, 



C+ ||f^m||r ^ Ir{Urn) ^{,Pm^ Ir{Um) 
2-n 



[f{\U^{z)\'')\Um{z)? - F{\Urn{z)?)\ dz 

> r{fi-l)F{\U^{z)\^)dz-27rr (11) 
Jo 



> (/^ - 1)C2 / 

Jo 



Um(z)\^^'+^^ - 1] dz-2TTf 



> - 1)C2 (ll^/mlliJliV) - C6 



for a constant cq > 0. 

By following the same arguments, we derive 

l^\\Ul,^fr < \\P^VIm{U^)\\ ■ ||f/l,m||r + 2 / f {\Umiz)\^)\U^{z)\\Ui,^{z)\dz 







< \\Ul,m\\r + 2C7 r {\UUz)\^'+^ + l) \Ui,Uz)\dz 

Jo 

< ||t^l,m||r + 2C7 II It^ml^'*'^"^ + l|| ^(s+i) || f/i m || L2(=+i) 

M 2s+l 

< ||^l,m||r + 2C7 (||C/'m||?''2't+i) + l) ||f^l,m||r 



and hence 



llC^i ,m\\r ^ Cg (||CAti||^2(s+i) ~I" l) • 

Similarly we obtain 

l|f^2 ,m||r ^ Cg ( II f^n 11^,2^+1) + l) 

and consequently by (fTT|) . we obtain 



\\Um\\r < 2C8 {WU^WltU + 1) < C9 (||f/m||^+^' + l) 

for positive constants cy, cg and cg. Thus {t/m}meN C Xr is bounded. Since 
and iCr^r is compact, there is a convergent subsequence of {?7m}meN ii^ 

. Sum- 
marizing, (PS)*-condition is verified for J^, i.e. condition (/i) of Theorem 12.31 is 
verified. 
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Next, let U G X„. Then using Ui £ spanjei, 62, • • • , e^}, we derive 

< \ - - C2 £^ - 1) dz 

< ^ (rill? - ^211.^) - cn (ll^ilir^^ + ||f/2|ir'0 + 
<^||C/i||?(l-ci2||C/ir)-^||C/2||? + cio 

for positive constants cio, cn and ci2- Now it is clear that Ir{U) —00 as \\U\\r 
00, i.e. condition (Is) of Theorem 12.31 is verified. 

Summarizing, under assumptions (HI) and (jlOp . all conditions (/i)-(/4) of The- 
orem [231 are verified for J^- Hence there is a nonzero critical point Ur G Xr of I,., 
which we already know to be a C^-smooth solution of ([6]) satisfying ([5]). Note (llOp 
certainly holds for any |z/| > ^ and r G (0, 1). Hence the proof of the second part 
of Theorem 11.11 is finished. To prove the first part, it enough to observe that the 
set 

{$(r + /c) I r G Qn (0,1), keZ} 

is countable, and thus for almost each G R \ {0} and any r G Q Pi (0, 1), condition 
(fTO|) holds. 



2.3. Remarks 

Remark 2 : When r is rational in Theorem 1 1 . 1 1 then we get periodic U{z) with 
arbitrarily large minimal periods. If r is irrational then clearly U{z) = e~^^^V [z) 
for a T-periodic V{z) = U{z)e~~'^^^ . So U{z) is quasi periodic and its orbit in 
C ~ IR^ is dense either in a compact annulus or in a compact disc. But |?7(-z)| is 
T-periodic in the both cases. 

Remark 3 : Changing t <-> —t, we can also handle DNLS 

-lUn = '^aj/\jUn + f{\Un\'^)Un-, nGZ (12) 

under (HI) and ([TO]) becomes 

1/ / $(r + A;)VA; G Z and i/G(0,^). (13) 
Remark 4 : Assume that U G is a weak solution of ^ , then 

\u{z)\ < \Uk\ < lY.w(^+^)\lY.^'+^y^ 

— \ — cosec vrr \\U \\y^. 
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Let R := maxx^iR^ x^{x). Then 

W\\\U'\\l^ = < W^rUU- + \\f{\U\^U)\\L- 



<R\\U\\l'^+ci\\\U\^'+^ + \U\\ 



< \^Rrl + cirl + ci^^cosec''^'TTr\\U\\f^j \\U\\y^ 
So if [/ 7^ then we obtain 

\u\ < Rrj + cirf + ci -^00860^"*+^ vrr 

i.e. 



TT' 



s+1 



U\\y >V2' 



' lu] — Rrf — cirf 
ci7r''+^cosec^''+^ vrr 



for 



\iy\ > Rr\ + cir\. 



Hence ||?7||y^ — > oo as — > oo for a possible nonzero solution U ^Yr oi ^ 



2.4. Examples 

We first note 



-y^Oj(l - cosxj) = - 

X ^ X 



jeN jew 



(14) 



Now we turn the the following concrete examples. 

Example 2.6 First we suppose that aj is decaying rapidly to 0. Let aj = jj. Then 



1 = e' 



cosx+ismx 



jew ■ 



gcos X g-^^ X + z sin sin xl — 1. 



So by (fn|) we derive 



— - e™*^^ cos sin X + 1 



- [e - e™*"^ cos sin xl 



By Remark [H $ € C^(IR, R) with the graph on [-47r,47r]: 
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3 

2 
1 




A numerical solution shows that $ has a maximum R = $(xo) = 3.15177 at 
xq = 1.03665. 

Example 2.7 Now we suppose that aj is decaying exponentially to 0. Let aj = 
, hence we have the discrete Kac-Baker interaction kernel [13, [3] • Then 



e^i' = e(^*-^)^' = - 



^Xl—l 



cos X + I sin X e cos x — 1+ ei sin x 



e — cos x — I sin x e"^ + 1 — 2 e cos x 



So by we derive 



E 



e cos X — 1 
e^ + 1 — 2 e cos x 



2e(e + 1)(1 - cosx) 
( e — l)x( e^ + 1 — 2 e cos x) 



By Remark [U $ G C^(IR, R) with the graph on [-47r,47r]: 





A numerical solution shows that $ has a maximum R = $(xo) = 0.992045 at 
xo = 0.991541. 

Ex amp le 2.8 In this example, we suppose that aj is decaying polynomially to 
(cf. [301]), by considering several cases: 

1. Let a, = Then 



1 1 

jew •' 



cos X] 



x\ — 2tx 



\A 

2-K 



2Ax 



271 - \x\ + 2tt 



2tt 



Here [•] is the integer part function. By Remark [H <I> € C^(IR,IR) with the graph 
on [— 47r, Att]: 
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$ has a maximum R = ^{xq) = ^ = 1.53117 at xq = 2tt/3 = 2.0944. Similar 
results hold for aj = with /5 > 3 by Remark [Tl 



2. Let aj = js- So we consider the dipole-dipole interaction (cf. [1, 13. [25l. [3o|). 
By RemarklH ^> G C{R, IR) with the graph on [-An, An]: 




$ has a maximum R = ^{xq) = 1.68311 at xq = 1.76076. Next we know that |4l(] 



- cos xj = — In 



2 sin ■ 



< X < 27r. 



Then 



sm 



1 In 


s 


2 sm — 




2 



Using x/2 < sinx < x for x > small, we derive 

X X 
X — xlnx = — / lnsds< sin xj < — / In — = x — x In — . 



jeN 
By L'Hopital's rule, we obtain 



_ $(x) _ 4E,,y^sm2xj 



lim — '^-^ = lim 

a;^0+ X x—*0+ 



X^ 



,. 2^ .g^ Jjsin2xj 

hm = +00. 

x^0+ X 



Hence $ has no derivative at xq = 0. 

Next, let aj = for 2 < /3 < 3. By Remark [Tl <I> is still continuous. Since 
^>(0) = and 

$(x) ^ ^. 4^.g^^sin2xj 
hm > hm ^ = +oo, 

3:^0+ X x^0+ X^ 
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<I>(x) is continuous but not C^-smooth on IR. 



3. Let aj = -k. Then 

3 J-! 



15 



x\ — 2tt 



2x 



2tt - \x\ + 27r 



\x\ 
2^ 



By Remarkm ^> G C{R \ {0}, R) with the graph on [-47r,47r]: 




$ is discontinuous at xq = where it has a supremum R = tt. 
4. Let Gj = for 1 < /? < 2. For (3 = 7/4, <I> has the graph on [-Air, Air]: 




Hence $ is discontinuous at xq = with hm2._jo+ ^{x) = +oo. We show that this 
holds for any 1 < (3 < 2. First suppose 3/2 < (3 < 2. Then the series 



sm jx 



converges uniformly on any [e, 27r — e] for < e < tt. But since ^jgD^ j2w-i) < oo, 
so T G C L^. On the other hand, we know 4]| that 



T(x) := r(2 - P) cos !^i/L_il . a;/3-2 + 0(1) 



on (0, vr]. Hence 



jew 

on [0,7r]. Consequently, we obtain 



^, , 2r(2-/3) ^..x 
$(x) = ^_ cos ^^^^ — ^ • x'^ 2 + 0(1) 
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on (0,7r], which imphes liuix^o^ ^{x) = +00 for any 3/2 < (3 < 2. Finahy, if 
1<P< 3/2, then 

, , , 2 v-^ 1 — cos jx 8 ^ / 1 \ 3-K 1 ^ , , 

as x — > 0+. Hence, Hm2.^o+ ^(x) = +00 for any 1 < /3 < 2. 

Summarizing, we have the fohowing result. 
Lemma 2.9: Let aj = j^^ for 1 < /?. Then 

(i) ^ e Ci(IR, R) /or /3 > 3, and 7^^> = [-R, R] for some R<oo. 

(ii) $ G C(IR, R) and $ ^ C^(R,IR) for 2 < p < 3, and = /or some 
R< 00. 

(Hi) $ e C(R \{0},IR) and $ ^ C(R,IR) for 13 = 2, andU^ = (-7r,7r). 

(iv) $ G C(R \ {0}, IR) and $ ^ C(R, IR) /or 1< /3 < 2, and = (-00, +00). 

Remark 5 : We see that if the interaction is strong, so the case (iv) of Lemma 
12.91 holds, then there are continuum many quasi periodic traveling wave solu- 
tions U{z) of Theorem 11.11 for any i/ 7^ 0, T > and r S (0,1) such that 
r ^ {z — [z] I z E ^<^~^(— I/)}, with ||C/||y^ ^ 00 as |z^| — > 00 by Remark|4l On the 
other hand, if the interaction is weak, then we can show in addition quasi periodic 
traveling waves with speeds in intervals (—00, —R) and {R, oo) for any T > and 
r G (0,1). 

Remark 6 : For the reader convenience, we present the above graphs of function 
$ to visualize their quantitative and qualitative changes according to different 
choices of values of sequences {ajjjgz in ©J and hence with different consequences 
from Theorems 11.11 and 11.21 for the existence and bifurcations of quasi periodic 
traveling wave solutions of ([3]). Moreover, these graphs can be compared with 
similar ones for traveling waves for higher dimensional DNLS in Section S] and for 
traveling waves with frequencies in Section El Finally these examples are motivated 
by applications mentioned in the corresponding references. 



3. Bifurcation of Traveling Wave Solutions 

In this section we proceed with the study of ([6]) when nonresonance of Theorem 
11.11 fails, i.e. r G {fi,f2, • • • ,fm}. We scale in ([6]) the velocity by i^/(l + A) to 
get equation 

-uiU'{z) = {l + \) {Y,ajd,U{z) + f{\U{zt)U{z) \ , (15) 
View / 

where A is a small parameter, i.e. Un{t) = U (^n — j^t^ is a solution of ([3|). We 

are interested in the existence of quasi periodic solutions U{z) of (fT5]) stated in 
Theorem 11.21 

3.1. Preliminaries 

In this subsection we recall some results from critical point theory of [29]. Let H 
be a Hilbert space with a scalar product (•,•) and the corresponding norm || • ||. 
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Let : 5^ — > L{H) be an isometric representation of the unit circle over H, 
i.e. the fohowing properties are satisfied 

(R) 6(0) = / - the identity, 9(^1 + 62) = 6(01)0(^2) for any 01,92 G S\ {e,h) 
@{6)h is continuous, and ||0(6')/i|| = \\h\\ for any 6 £ and h £ H. 

We set 

Fix(5^) ■={heH\ Q{e)h = h\f9 G G} . 

We consider Ji, J2 G C^iH, R) such that 
(HI) J2(0) = and VJi(O) = VJ2(0) = 0. 

(H2) Hess Ji(0) is a Fredholm operator, i.e. dim Hess Ji(0) < cxd, T^Hess Ji(0) is closed 

and codimT^Hess Ji(0) < 00. 
(H3) dim ker Hess Ji (0) > 2 and Hess J2(0) is positive definite on ker Hess Ji(0). 
(H4) Ji and J2 are S'^-invariant, i.e. Ji,2{®{G)h) = Q{6)Ji,2{h) for any G and 

h£H. 

(H5) ker Hess Ji(0) nFix(5i) = {0}. 

Now we can state the following [2^ . Theorem 6.7]. 

Theorem 3.1: Under the above assumptions (H1)-(H5), for each sufficiently 
small e > 0, equation 

VJi(/i) + AVJ2(/i) = (16) 
has at leat | dim ker Hess Ji (0) of S^- orbit solutions 

{{\k{e),Q{e))hk{e) | G 5^} , A; = 1,2,--- , ^ dim ker Hess Ji (0) 

such that J2{hk{e)) = e and hk{e) 0, Afc(e) ^0 as e — > 0. Clearly hk{e) 7^ 0. 

Remark 1 : When Hess J2 (0) is negative definite on ker Hess Ji (0) , then Theorem 
13. II holds for e < small. 

Remark 2 : By (H4), ker Hess Ji(0) is invariant with respect to 0. Using (H5), 
dim ker Hess Ji (0) is even. 

Now assume H = H^(BH^ be an orthogonal and 0-invariant decomposition with 
the corresponding orthogonal projections P± : H H±. Then Q(9)P± = P±@{9) 
for any G 0. Let us consider an equation 

C(/+ - I^)h + (1 + A)(0 + VJ^ih)) = 0, (17) 

where C 7^ is a constant, A is a small parameter, I± : H± H± are the identities. 
We suppose 

(A) }C:H^His compact self-adjoint and T G C^iH, R) with J^(0) = 0, V:r(0) = 0, 
Hess .7^(0) = 0, and /C, are 5^-invariant. Moreover, ]CH± C H±. 

Then 

j,{h) = ^-{\\p+hf - WP^hf) + i(o, h) + T{h), 

J2{h) = ^ilCh,h)+J^{h). 
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Hence 

Ji(0) = J2(0) = 0, VJi(O) = VJ2(0) = 0, 
HessJi(O) = C(/+ - /-) +/C, HessJ2(0) = /C. 

So assumptions (HI), (H2) and (H4) are satisfied. Since P±IC = ICP±, equation 

Hess Ji(0)/i = C{I+ - I-)h + lCh = 

splits into 

Consequently, supposing either 

(B+) ker(C/ + JC)nH+ = {0}, dimker(CI - JC) n > 2 and ker(Ci" - /C) n n 
Fix(5i) = {0} 

or 

(B_) ker(C/ - /C) n = {0}, dimker(CI + JC) n H+ > 2 and ker(Ci" + IC) n H+ n 
Fix{S^) = {0} 

we get either 

ker Hess Ji (0) = ker(C/ - /C) n F_ 

or 

ker Hess Ji (0) = ker(C/ + fC) n H+ 

and so (H5) holds as well. Finally, we derive 

Hess J2 (0) I ker Hess Ji (0) = ±(1 

and thus (H3) is also verified (cf. Remark [1]). Summarizing, Theorem 13.11 and 
Remark [1] is applicable to (fTTj) : 

Corollary 3.2: Under assumptions (A) and (B±), for each sufficiently small 
e 7^ 0, ±e(" > 0, equation (fTTl) /las at leat ^ dim ker (^J =F /C) H ffzp of -orbit 
solutions 

{{\k{e)MO))hk{e) \e&S^], k = l,2,--- , ^ dimker(CI T /C) n 

such that ^(/C/ifc(e), /ifc(e)) + J^(/ifc(e)) = e and /ifc(e) 0, Afc(e) ^ as e ^ 0. 
Clearly /ifc(e) / 0. 

Remark 3 : If Fix(S'i) = {0} then (B+) holds if 

(i) — C ^ (t(/C///+), C G (y{K./H^) and C has a multiplicity at least 2, 
while (B_) holds if 

(ii) — C G (t(/C///+), C ^ (t{K./H^) and — C has a multiplicity at least 2, 
respectively. 
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3.2. Proof of Theorem [TTH 

We again assume for simplicity T = 2tx. So let r = fi S (0, 1) and the equation 

— ly = (J) (fi + k) 

has solutions ki,k2, - ■ ■ , kmi G 2 which are either all nonnegative, or all negative. 
Next ([I5|) has the form (cf. ([9])) 

2{iyl+ - !//_) - (1 + A) (krCrU + ^r(f/)) = (18) 

and 

IC = -KrCr, T{u) = - F{\U{z)\'^)dz. 

Jo 

Isometric representation is naturally given as 

@{e)u{z) := u{z + e), 

i.e. 

\kez / kez 

Note Fix(5^) = {0}. It is easy to verify (R) for 6. By results of Section [2l we get 
both ICH± C H± and assumption (A) holds, and moreover 

a {}C/H±) = {±2<^{n + k)\k£Z±}. 

Note Z+ = {0} U IhJ and Z_ = -\s4. Hence (i) of Remark [3] is satisfied if 

-V ^ {$(fi + k) \ k £ Z+} , -u E {^>(fi + k) \ k e Z_} , 

while (ii) if 

-1/ G {$(fi + A:) I /c G Z+} , -u ^ {$(fi + /c) I /c G Z_} . 

But these are precisely assumptions of Theorem 11.21 So its proof is complete by 
Corollarv 13.21 and Remark [3l 
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4. Traveling Waves for Higher Dimensional DNLS 

In this section, we first show how to extend previous results for 2-dimensional 
DNLS (2D DNLS) [13, [HI, HI] of forms 



(19) 



where «n,m G C, := \ {(0,0)}, AijUn,m ■= Un+i,m+j + Un-i,m-j " 2tt„,m are 
2-dimensionaI discrete Laplacians, / satisfies (HI) and Uij = a^i-j along with 
|ai jl < oo and all Ui j are not zero. 

Again, (jl9p conserves two dynamical invariants 



(n,m)S 



|ttn,m|^ — the norm, 

(n,m)eZ2 

^ ^ Ojj "Urijml ~l~ -^(|'Wn,m| ) 



the energy. 



We look for traveling wave solutions of (jl9p of the form 



Un,m{t) = U{n COS 6 + 171 sin — vt) 



(20) 



with a direction (cos 6*, sin 0) [2l|]. Hence we are interested in the equation 



mU'iz)= a^J^i,JUiz) + fi\U{z)\^)Uiz). 



(21) 



where z = n cos 6 + m sin — i^t, 7^ and 



dijU{z) := U{z + i cos 6* + j sin 6*) + U{z - icos9 - j sinO) - 2U{z). 



We see that (|21|) has a very similar form like ([U]). So we can directly repeat the 
above arguments, where now instead of ^(x) we get 



$e(x) 



2 x(z COS + j sin ( 



Set Rg := sup|]^<I>0. Summarizing, Theorems 11.11 and 11.21 have the following analo- 
gies: 

Theorem 4.1: Let (HI) hold and T > 0, 9 e [0, 27r). Then for almost each 
G IR \ {0} and any rational r G Q Pi (0, 1), there is a nonzero periodic traveling 
wave solution (I20p of (119p with U G C"'^(IR,C) satisfying ([5|). Moreover, for any 
v £ R \ {0}, there is at most a finite number of fi^g, f2^e, • ■ ■ , ^m£,,e G (0, 1) such 



September 10, 2009 



17:14 Applicable Analysis revisedlatticeFeckanRothos 



Traveling Waves of DNLS 



21 



that equation 



has a solution k ^ Then for any r G (0,1) \ {^i,6») "^2,6) " " " i^me,d} there is a 
nonzero quasi periodic traveling wave solution (I20p of (I19p with the above proper- 
ties. In particular, for any > Rg and r £ (0, 1), there is such a nonzero quasi 
periodic traveling wave solution. 

Theorem 4.2 : Suppose f G C^iR+, R) with /(O) = 0. // there are G (0, 1), 
T > 0, G [0, 27r) and v G Tl^e \ {0} such that all integer number solutions 
ki,k2, - ■ ■ , krrii a of equation 



(2-K 



are either nonnegative or negative, and mifi > 0. Then for any e > small there 
are mi^g branches of nonzero quasi periodic traveling wave solutions (j20p of (jl9p 
with Uj^e £ C^{R, C), j = 1,2, ■ ■ ■ , mifi, and nonzero velocity satisfying Uj^s{z + 
T) = e^'^'"^*C/(z)j^£, Vz G IR along with Vfr ^ v and Uj^e ^ uniformly on R as 



Example 4.3 We consider the discrete 2D Kac-Baker interaction kernel ajj- 
e-Nhlil for {i,j) G Zg. Then ^ e'l^l-l^'l 



4e 



- 

(e^ 



and 



1 



(e — 1)2 (1 + e^ — 2 e cos(x cos 0))(1 + e^ — 2ecos(2;sin6')) 



A numerical evaluation shows that function {x,6) ^9{x) has a maximum R = 
9.75047 at xq = 1.08205 and 6*0 = 0.785398. To justify this theoretically, we take 
a = xcosO and b = xsinO to transform ^e{x) into 



Ha,b) 



{e + lf (e^-l)^ 

(e-l)2 (1+ e2 - 2ecosa)(l + e2-2ecos6) 



Va2 + 62 



Note $(a,6) = <I>(iba,ib6) = $(5, a). A numerical evaluation shows that function 
^s{a,b) has a maximum R = 9.75047 at oq = feo = 0.765123 which correspond to 
xo and 6^0. On the other hand, if a2 + ^2 > 4 then $(a, b) < 2^^^ = 9.36539 < 



9.75047, so $(a, 6) achieves its maximum in the disc D2 
solving the system ■^^{a,b) = ^$(a,6) 



{a2 + ^2 < 4}. Next, 
at the max- 



we derive - -f- 

ao bo 

imum point (ao,6o) G D2, oq > 0, bo > 0. But the function ^^^^ is decreasing on 
[0,2], so oq = bo, and thus 
for function 



7r/4. An elementary but awkward calculus shows 



(e + 1) 



(e2 



1) 



1)2 



1 + e2 — 2 e cos 



{.4 



4 

X 



with the graph on [—20, 20]: 
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that Xq £ (0, 2) is the only root of ^'^^^(xq) = on (0, 2), and then R = ^^^^^{xq). 

So R is computed also analytically in this case. 

Summarizing, Theorems 14.11 and 14.21 can be applied in this case for any suit- 
able nonzero z^, and resonant traveling waves with maximum velocities which are 
achieved in the diagonal directions it^o = ±vr/4. 

Finally, it is now clear how to proceed to 3D DNLS or even to higher dimensional 
DNLS, so we omit further details. 



5. Traveling Waves with Frequencies 

We could consider more general traveling wave solutions than above of forms 

ujt) = Uin-vt) e*"^*, 

(22) 

Un,m{t) = U{ncos6 + msinO - ut) e*"^* 



with velocity ^ and frequency lo ^ (see [33] )• Then, there is a dispersion 
relation between the velocity and frequency to as follows. Inserting (122p into ([3]) 
and (fT9jl . respectively, we are interested in equations 

-mU'iz) = ajdjU{z) + ujU{z) + f{\U{z)\^)U{z), 

(23) 

-mU'{z) = ai,,dijU{z) + ujU{z) + f{\U{z)\^)U{z), 

respectively. We see that ([6]), (f2T]) and (p3]) are very similar, so we can repeat the 
above arguments to (j23p when instead of <I>(x) and ^e{x) now we have 

^'(x, oj) := -^{x) - - , $e(x, uj) := <^0{x) - - , (24) 

X X 

respectively. Consequently, we have analogies of Theorems 11.11 11.21 14.11 and 14.21 to 
(I23p but we do not state them since they are obvious. 

Example 5.1 We consider the discrete Kac-Baker interaction kernel from Exam- 
ple O Then 

^, , 2e(e -I- 1)(1 — cosx) uj 

$(X, LO) = . , „ . 

(e — l)x(e^ -|- 1 — 2ecosx) x 

To be more concrete, we first take w = 1, and then <I>(x, 1) has the graph on 
[-47r,47r]: 
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04 - 




-0.4 - 



with lim^^Oi ^{x,\) = =fcxd. A numerical evaluation shows that function ^>(3;, 1) 
has a maximum R = 0.282071 on (0,oo) at xq = 1.9905. Consequently, the analogy 
of Theorem 11.11 can be applied now to any v ^ {) while the analogy of Theorem 
Ocan be applied for almost any i/ G R \ [-0.282071,0.282071], while for nonzero 
V e [-0.282071,0.282071] could be problematic in general. 

On the other hand for lo = —1, —1) has the graph on [— 47r,47r]: 




with lim^^^Oi ^{x, —1) = itcxD. Consequently, the analogy of Theorem 1 1 . 1 1 can again 
be applied now to any 7^ while the analogy of Theorem 1 1 . 2 1 can now be applied 
for almost any v ^ Oi course now we have totally different situations than 
in Example 12.71 for traveling waves without frequencies by comparing the above 
graphs with that one in Example 1 2. 7[ 
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